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ABSTRACT 


This  report  consists  of  tvo  separate  articles  in  which  microwave  analog 
measurements  of  scattering  by  non- isotropic  spheres,  stacked  7-cylinder 
rough  particles,  and  2n(n-l,2,3)  sobers  are  presented.  Emphasis  is 
placed^on  extinction  measurements  te-Qo) scattering)  but  some  side  scattering 
data  (ji00<e<lUog)  is  presented  in  the  first  paper  which  is  devoted  solely 
to  scattering  By  multiple  spheres.  . 
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FOREWORD 


In  1958,  Dr.  J.  Majro  Greenberg  (now  at  the  University  of  Leiden, 
Netherlands)  established  a  microwave  analog  facility  to  investigate  the 
extinction  of  light  by  small,  irregular  particles.  Joined  by  Dr.  Ru  T. 
Wang  in  i960  and  later  by  Dr.  Donald  Schuerman  in  1972  (who  now  directs 
the  facility) ,  the  facility  has  fostered  numerous  investigations  of  both 
the  extinction  and  angular  scattering  of  radiation  by  particles  whose 
size  is  of  the  order  of  the  incident  wavelength.  Many  of  these  studies 
have  gone  unpublished.  Under  this  contract,  we  have  selected  a  fairly 
large  number  of  intriguing  scattering  results,  put  them  in  a  form  appro¬ 
priate  for  publication,  and  tried,  where  possible,  to  provide  physical 
interpretations.  Our  desire  is  to  disseminate  the  results  of  these 
unique  measurements.  This  report  is  part  of  that  endeavor.  It  consists 
solely  of  two  articles.  The  first,  entitled  "Extinction  Signatures  of 
Kon-Spherical/Non-Isotropic  Particles"  by  Ru  T.  Wang,  appeared  in  the 
book  LIGHT  SCATTERING  BY  IRREGULARLY  SHAPED  PARTICLES  (Edited  by  D.  W. 
Schuerman,  Plenum  Publishing  Corp. ,  1980).  The  second,  entitled 
"Scattering  by  Multiple  Spheres"  by  R.  T.  Wang  and  D.  W.  Schuerman,  will 
be  presented  at  the  1980  CHEMICAL  SYSTEMS  LABORATORY  SCIENTIFIC  CONFERENCE 
ON  OBSCURATION  AND  AEROSOL  RESEARCH  and  will  hopefully  be  published  as  a 
part  of  those  proceedings.  An  abbreviated  or  summary  version  of  these 
articles  will  soon  be  submitted  for  publication  in  the  journal  APPLIED 
OPTICS  (reprints  of  which  will  not  be  available  for  many  months).  That 
summary  article  will  refer  interested  readers  to  the  original  articles  and 
to  this  report. 

Don  Schuerman 

Albany,  NY 

July,  1980 
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SCATTERING  BY  MULTIPLE  SPHERES 


R.T.  Vang  and  D.W.  Schuerman 
Space  Astronomy  Laboratory 
State  University  of  New  York  at  Albany 
Albany,  NY  12203 


ABSTRACT 

Microwave  analog  measurements  of  scattering  by  ensembles  of  2n(n=l,2,3) 
interacting  spheres  are  presented.  Emphasis  is  placed  on  the  amplitude  and 
phase  measurements  of  the  forward  (8=0°)  scattered  wave  as  a  function  of  the 
mutual  separation  of  the  spheres  and  the  orientation  of  the  ensembles  with 
respect  to  the  incident  beam.  The  results  are  displayed  as  calibrated  P,Q 
plots;  i.e.,  cartesian  displays  of  the  complex  amplitude.  A  short  historical 
sketch,  a  description  of  the  experiment,  and  phenomenological  explanations  of 
the  results  are  also  given.  Side  scattering  (0=90°)  and  angular  distri¬ 
bution  (J»0°sesl40°)  measurements  indicate  the  existence  of  specular  scattering 


1.  INTRODUCTION 


Small  particles  scatter  light  quite  differently  depending  on  vhether 
they  are  mutually  veil  separated  or  in  close  proximity.  This  is  quite 
obvious  if  ve  consider  the  tvo  extreme  cases  of  2  spheres  either  far  apart 
or  in  contact  (a  single  particle).  The  earliest  theoretical  investigation 
of  dependently  scattering  particles  seems  to  be  that  by  Trinks  (1935)  who 
formulated  the  solution  of  Maxvell's  equations  for  a  2-sphere  problem  of 
rather  small  particle  size.  Despite  the  mathematical  complexity,  theoretical 
work  continued  (see  Germogenova,  1963;  Liang  and  Lo,  1967;  Levine  and  Olaofe, 
1968;  Rozenberg,  1971;  Bruning  and  Lo,  1969,  1971).  In  these  rigorous 
solution  approaches,  one  still  sums  the  scattering  by  individual  particles 
to  evaluate  the  total  scattering  from  the  array.  To  compute  the  scattering 
from  each  sphere,  taking  into  account  the  near  field  effect  of  the  neigh¬ 
boring  particle  requires  the  rather  complex  procedure  of  relocating  the 
spherical  vector  wave  functions  from  one  particle  origin  to  the  other  to 
form  a  proper  boundary-value  problem.  This  procedure  results  in  extremely 
complicated  expressions.  Although  the  advent  of  modern  computers  and  the 
development  of  efficient  algorithms  make  the  numerical  evaluation  more 
accessible  (Bruning  and  Lo,  1971),  it  is  still  difficult  to  obtain  a  clear 
physical  picture  of  dependent  scattering.  We  have  also  to  mention  the  very 
ingenious,  yet  complex  theoretical  vorks  of  Tversky  (1967),  Waterman  and 
Truell  (1961)  and  the  less  formidable  but  still  difficult  vorks  of  Hongo 
(1978),  Borghese  et  al.  (1979)  and  Kattavar  and  Humphreys  (1980).  Like 
those  in  the  single-particle  cases  (van  de  Hulst,  1957;  Kerker,  1969),  the 
problem  of  multiple-particle  scattering  has  many  stumbling  blocks  if  pur¬ 
sued  by  theoretical  approaches  alone.  Indeed,  even  in  a  static  field  the 
problem  of  finding  the  induced  dipole  moment  on  2  spheres  as  a  function 
of  mutual  separation  is  already  quite  involved  (Goyette  and  Navon,  1976). 

Experimental  investigations  on  multiple-particle  scattering,  although 
scarce,  have  also  been  sporadically  reported  during  the  past  tvo  decades. 
Perhaps  such  a  study  is  next  to  impossible  in  the  optical  region  (Woodvard, 
196U),  and  the  few  experiments  conducted  in  the  microvave  region  vere  per¬ 
formed  mostly  for  the  backscattering  (0=l8O°).  Angelakos  and  Kumagai 
(196U)  obtained  backseat ter  results  for  multiple  conducting  spheres  and 
compared  them  vith  the  predictions  of  geometrical  optics.  To  verify  their 
theoretical  results,  Bruning  and  Lo.  (1969,  1971)  performed  such  experiments 
for  both  conducting  and  penetrable  multiple  spheres.  The  only  exception 
in  vhich  an  extinction  (0«O°)  experiment  vas  performed  is  by  Beard,  et  al. 
(1967  and  references  cited  therein)  vho  investigated  multiple  spheres  in 
random  motion,  recorded  the  quadrature  phase  components  of  such  aggregates, 
and  subsequently  compared  the  results  vith  the  statistical  theory  of  Havley 
et  al.  (1967).  Disregarding  the  technical  difficulties  of  using  millimeter 
vaves  for  the  0*0°  scattering  research,  the  random  particle  motion  alone 
may  completely  obscure  the  detailed  picture  of  dependent  scattering. 

Encouraged  by  the  initial  success  of  0=0°  microwave  measurements 
(Lind  et  al. ,  1965)  for  non-spherical  particles  vhich  permitted  us  to  ob¬ 
serve  detailed  phenomena  vith  precisely  known  target  parameters  and  orien¬ 
tations,  ve  extended  the  analog  method  to  cover  the  dependent  scattering 
studies  by  replacing  a  non-spherical  particle  vith  an  array  of  spheres. 

Some  0*90°  and  l»O°<0<ll*O°  scattering  measurements  vere  also  made  to  explore 


the  elegant  mathematical  symmetry  properties  and  specular  scattering  phenomena 
exhibited  by  such  an  array.  Our  motivation  was  also  prompted  by  the  afore¬ 
mentioned  scarcity  of  relevant  data  and  partly  by  the  encouragement  of  the 
late  Prof.  P.  Debye  who  recognized  the  significance  of  such  experiments. 

The  data  presented  in  this  report  were  accumulated  during  the  period  1968- 
1970  in  two  different  buildings  at  Rensselaer  Polytechnic  Institute,  Troy, 

New  York.  \ 

2.  ON  PERFORMING  THE  EXPERIMENTS 

In  former  theses,  reports,  and  published  works  (Lind,  et  al. ,  1965;  >- 

Lind,  1966;  Wang,  1968;  Wang  et  al.,  1977;  Wang  and  Greenberg,  1978)  we 
have  already  described  the  experimental  method  in  great  detail.  In  this 
report  we  only  supplement  these  descriptions  with  topics  not  previously 
covered  or  those  important  in  performing  multiple-sphere  scattering  ex¬ 
periments. 

2.1  Target  Preparation  and  Refractive  Index  Determination 

The  preparation  of  scatterers  of  accurately  known  size,  shape  and  re¬ 
fractive  index  is  an  essential  step  in  performing  a  microwave  scattering 
experiment.  The  scattering  targets  are  manufactured  by  either  molding  or 
machining  commercially  available  plastic  materials.  Since  the  refractive 
index  of  such  a  material  depends  only  on  its  density  under  normal  laboratory 
conditions,  the  molding  technique  provides  an  adequate  means  of  obtaining 
a  scatterer  of  the  desired  refractive  index.  The  expandable  polystyrene 
supplied  by  the  Sinclair-Koppers  Co.  under  the  trade  name  Dylite  F-40  was 
found  to  possess  stability,  mechanical  strength,  low  humidity-absorption 
and  conductivity.  The  commercial  supply  comes  in  small  beads  about  0.5  mm 
in  diameter,  each  with  impregnated  volatile  material  to  help  expansion 
when  heated.  A  proper  amount  of  such  beads  will  expand  and  fuse  to  each 
other  when  heated  in  an  enclosed  cavity,  transforming  them  into  a  strong, 
smooth-skinned  foam  filling  the  mold  cavity.  Three  differently  sized  spherical 
cavities  (3.2  cm,  3.8  cm,  and  4.7  cm  diameters)  in  separate  stainless  steel 
blocks  were  used  for  molding  to  facilitate  the  mold- release  and  to  allow 
heating  by  steam.  Preweighted  amount  of  beads,  enough  for  a  near-tight- 
fill  was  poured  into  each  mold  and  was  steam  heated  for  about  5  hours.  Slow 
expansion  in  the  molding  process  was  preferred  to  insure  the  homogeneity 
of  the  product  target  medium.  At  least  9  identical  particles  for  each  size 
were  thus  fabricated  to  conduct  the  multiple-particle  experiment. 

The  refractive  index  of  the  target  was  determined  through  the  di¬ 
electric  constant  measurement  of  rectangular  wave  guide  samples  prepared 
from  the  same  Dylite  F-40  molded  in  short  pieces  of  the  waveguide.  The 
density  of  such  samples  spanned  that  of  the  target  sample  so  that  the  re¬ 
fractive  index  of  the  latter  could  be  linearly  interpolated  against  the 
density  through  the  measured  values  of  the  former.  The  measurement  tech¬ 
nique  employs  the  classical  standing  wave  method  originated  by  Roberts  and 
von  Hippel  (1946)  and  its  development  in  this  laboratory  (see  also  Sucher, 

1963;  Westphal,  1954).  The  conceptual  simplicity,  accuracy,  and  versatility 
of  the  method  allows  the  following  condensed  description. 

A  standing  wave  is  set  up  inside  a  waveguide  by  an  incident  wave 
traveling  toward  an  impedance  discontinuity  (the  shorted  waveguide  terminal) 


it 


' 


and  a  reflected  wave  from  the  same  discontinuity.  The  positions  of  miniffia/ 
and  the  amplitude  of  the  standing  wave  depend  primarily  on  the 
geometry  of  the  guide  and  the  dielectric  property  of  the  medium  containing 
this  wave.  Hence,  by  comparing  the  standing  wave  patterns  with  and  without 
the  dielectric  sample  in  the  guide  we  can  deduce  the  complex  dielectric 
constant  e  *  e*-ie"  of  the  sample.  The  complex  refractive  index  m  «  m'-im" 
then  follows  Maxwell's  relation 


or  equivalently 


m»2-m"2  *  c* 


2m' m"  =  e" 


Figure  1  is  a  schematic  of  the  dielectric  measurement.  A  traveling  probe 
(HP  ,  fitted  to  a  slotted  section  of  a  rectangular  wave  guide  (HP  X809B), 

measures  the  standing  wave  profile.  The  position  of  the  probe  along  the  guide 
is  continuously  monitored  by  a  dial  gauge  (Ames  282M)  to  an  accuracy  of  ±0.001  cm. 
One  end  of  the  slotted  section  is  short-circuited  by  a  flat  brass  plate;  the 
other  end  is  connected  in  series  to  the  klystron  oscillator,  isolator,  frequency 
meter  and  a  precision  attenuator  (HP  X382A).  The  oscillator  is  repeller  modu¬ 
lated  to  produce  a  1000  HZ  square-wave  modulated  microwave  of  carrier  frequency 
fQ  =  9. !»17  GHZ,  the  same  frequency  used  in  the  scattering  set  up.  The  electric 
field  of  the  standing  wave  is  read  by  the  standing  wave  indicator  (HP  4l5B) 
connected  to  the  probe.  Let  Emi n  and  Emnr  be  the  minimum  and  maximum,  res¬ 
pectively,  of  such  a  field.  A  rectangular  dielectric  sample  of  length  1£ 
snuggly  fits  the  guide.  Let  x©  be  the  distance  between  the  sample  surface  and 
the  adjacent  E^n  position  in  the  empty  portion  of  the  guide.  The  other  sample 
face  touches  the  short-circuit  plate.  The  propagation  constants  within  the 
guide  are  denoted  by  Y;  and  Y2  Tor  the  empty  protion  and  the  sample-filled 
portion,  respectively.  Let  XQ  be  the  free-space  wavelength,  Ag  be  that  in¬ 
side  the  guide  of  width  a  and  height  b,  while  Ac  =  2a  is  the  cut-off  wave¬ 
length  of  the  same  guide.  Also  let 

k  <*  2ir/A  ,  k  *  2tt/A  ,  k  ■  2ir/A  .  (2) 

o  o’  g  g’  c  c 

If  the  empty  waveguide  is  sufficiently  loss-less  and  the  sample  is  nonmagnetic, 
an  analysis  of  the  transmission  of  the  TE01  mode  microwave  shows  (Roberts  and 
von  Hippel,  19^6;  Westphal,  195**): 
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STANDING  WAVE  INDICATOR 


and  therefore 


e 


k^ 


(6) 


Thus,  the  desired  e  is  found  by  (6)  in  terms  of  y%  which  in  turn  is  found  by 
(3)  through  the  measured  quantities  1e,  E^  n/Ernay  =  1/r  =  inverse  standing 
wave  ratio,  xQ,  and  kg.  In  practice,  the  slotted  line  is  first  short- 
circuited  without  the  dielectric  sample,  and  two  successive  probe  positions, 
and  Dpi,  for  Eg^  are  recorded.  From  these,  Ag  is  found  by 


(T) 


The  sample  is  then  inserted  into  the  guide  with  one  of  its  end  surfaces  squarely 
touching  the  short-circuit  plate.  A  new  Eg^  position,  D,  is  found  such  that 


D„>  D  >  D_ 


(8) 


Since  xQ  occurs  only  as  tan(kgxQ)  in  (l),  a  straightforward  consideration  yields: 


tan(k  x  )  =  -tan 
g  o 


yv1*1.1 


t: 9 ) 


which  replaces  the  measurement  of  xD  by  that  of  DR,  D  and  le.  The  measurement 
of  r  =  Emax/Emin  carried  out  ^y  means  of  the  calibrated  variable  attenuator 
(HP  X  382A)  from  which  AR(in  dB),  the  difference  in  attenuation  to  bring  Egp-n 
and  to  the  same  output  level,  is  obtained.  AR  is  directly  related  to  r 

by 


r  =  io.UR/2°.) 


(10) 


This  method  has  the  advantage  of  being  independent  of  the  detector  characteristics 
of  the  probe.  With  these  measured  parameters,  (3)  is  solved  numerically  for 
Yjl  by  the  Newton-Raphson  interation  method.  There  are  an  infinite  number  of 
possible  roots  for  Y2le  arising  from  the  multiple-valued  inverse  tangent 
function.  The  true  root  is  selected  from  the  two  sets  of  roots  obtained  from 
two  samples  of  the  same  material  but  with  different  le's.  Only  the  correct 
root  appears  in  both  sets.  Better  still,  if  an  approximate  value  of  msm'-im" 
is  known  and  m"  is  small,  a  good  starting  point  (C0,n0)  for  the  above  iteration 
in  the  complex  plane. 
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is  to  pick 
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V^)!-  &2 . 


(12) 


a  logical  choice  in  view  of  (5).  With  the  desired  Y21e  and  hence  Y2»  e  and  11 
follow  immediately  from  (6)  and  (l).  Finally,  the  measurement  error  correction 
on  the  clearance  gap  between  the  waveguide  dimension  b  and  the  sample  thick¬ 
ness  d  is  considered  based  on  the  reasoning  of  the  lines-of-force  distribution 
across  the  dimension  b.  The  corrected  value  ec  =  z'c  -  ie|t  is  given  by 


.  -  e'd 
ec  "  b-(b-d)e' 

(13) 

„  e"bd 

Ec  ~  ^b-(b-d)E*|jz 


2.2  Scattering  Quantities  and  Symmetry  Relations 

The  scattering  pattern  of  a  symmetric  particle  always  displays  symmetry 
with  respect  to  its  orientation  about  the  incident  beam.  A  careful  con¬ 
sideration  of  this  property  greatly  reduces  the  number  of  required  measure¬ 
ments.  We  have  already  discussed  this  for  the  0=0°  scattering  experiments 
(Wang,  1968;  Wang  and  Greenberg,  1978)  for  particles  of  rotational  symmetry, 
and  a  short  extension  to  other  scattering  angles  0  will  be  mentioned  for 
similar  particles.  A  pair  of  spheres  is  a  symmetric  "particle”  because  the 
axis  of  rotation  passes  through  the  spheres'  centers. 

All  single  particle  scattering  quantities  can  be  defined  through  the  use 
of  the  scattering  amplitude  matrix  §,  the  ^  elements  of  which  are  dimensionless 
complex  numbers.  This  is  explained  in  great  detail  in  van  de  Hulst's  text 
(1957),  and  we  adhere  to  his  notations  throughout  this  report. 

Figure  2  shows  the  scattering  geometry.  A  linearly  pglar^zed  incident 
wave  whose  electric,  magnetic  and  propagation  vectors  are  E0,  H0  and  k0,  re¬ 
spectively,  progagates  along  the  z-axis  of  a  laboratory  fixed  coordinate  frame 
(x,  y,  z)  whose  y-z  plane  is  chosen  horizontal.  This  plane  is  called  the 
"scattering  plane".  The  receiver  antenna  is  moved  through  this  plane,  at  a 
constant  distance  from  the  target  site,  to  observe  the  scattered  wave.  The 
polarizations  of  the  transmitting  and  receiving  antennas  were  kept  vertical 
(parallel  to  the  x-axis)  throughout  this  report.  The  pair  of  spheres,  witlj 
center-to-center  distance  s  and  whose  symmetry  axis  makes  an  angle  x  with  kD 
as  shown  in  Figure  2,  was  always  rotated  in  the  y-z  plane  (the  k-H  plane) 
during  an  experiment.  The  angle  x  is  called  the  "orientation  angle."  Th$ 
scattering  angle  0  is  the  angular  position  of  the  receiving  antenna  from  kQ 
direction,  and  at  which  the  observation  of  scattered  wave  is  made.  A 
straight  line  bisecting  the  supplement  of  0  in  the  y-z  plane  is  called 
"bisectrix",  and  a  plane  containing  this  bisectrix  and  orthogonal  to  the  y-z 
plane  is  called  the  "bisectric  plane". 
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Given  an  arbitrary  scatterer,  van  de  Hulst  (1957)  explained  that  there 
existed  four  related  particle  orientations  with  respect  to  the  incident  polari¬ 
zation  in  which  S  at  each  orientation  could  be  represented  by  the  same  set 
of  matrix  elements.  Si,  S2,  S3,  and  S4: 


Orientation  a  Orientation  b  Orientation  c  Orientation  d 

S  =/S2  S3\  Sv=/S2  "S,A  8  =(  S2  8.  =  /S2  S*»\ 

a  \S4  sj  b  \-s3  sj  C  V-S4  S J  d  \s3  sj 


Orientations  b,  c  and  d  are  achieved  from  the  orientation  a  through:  180° 
particle  rotation  about  the  bisectrix,  mirroring  the  particle  with  respect  to 
the  scattering  plane  and  mirroring  the  particle  with  respect  to  the  bisectric 
plane,  respectively.  For  a  rotationally  symmetric  particle  the  immediate  con¬ 
sequences  of  these  relations  are  the  following: 

(1)  If  the  rotation  axis  lies  in  the  scattering  plane,  the_ 
orientations  (a)  and  (c)  are  identical.  Equating  the  S 
matrix  elements  for  both  positions  we  have  S3=-S3  and 
si»=“S4,  which  is  true  if  and  only  if  S3=S4=0.  With  zero 
off-diagonal  elements  in  the  S  matrix,  there  is  no  cross- 
polarized  component  in  the  scattered  wave. 

(2)  As  the  symmetry  axis  is  rotated  in  the  scattering  plane, 
the  scattered  signal  varies  periodically  with  respect  to 
the  X  variation.  The  pattern  is  symmetric  about  the  bi¬ 
sectrix  with  period  n/2.  That  is,  the  symmetry  axis 
needs  only  be  swept  in  the  x  interval: 

f<x<f  +  f.  (HO 


This  follows  from  the  property  S 3=84=0  under  this  rotation. 

Outside  the  range  of  X ,  one  can  always  find  a  reciprocal 
orientation  position  inside  the  x  range  (equally  spaced 
from  the  bisectrix)  whose  corresponding  orientation  has 
the  same  S2  and  Sj. 

2.3  On  Suspending,  Orienting  and  Separating  Multiple  Spheres 

Since  the  objective  is  to  find  the  precise  scattering  signature  of  an 
ensemble  of  spheres  as  a  function  of  its  geometrical  configuration,  an  ideal 
suspension  and  orientation  mechanism  should  add  little  to  the  true  scattering 
from  the  ensemble  and  should  allow  the  quick,  accurate  and  reproducible 
positioning  of  the  ensemble  to  the  desired  orientation.  It  was  found  that 
a  skillful  employment  of  the  orientation  device  described  in  Lind's  thesis 
(1966)  with  only  a  few  modifications  would  meet  most  of  the  requirements. 


Two  small  nylon  eyelets  were  anchored  onto  a  pair  of  diametrical  poles 
of  each  sphere,  and  thin  nylon  strings  passing  through  each  eyelet  allowed 
the  assembling  of  multiple  spheres.  Appropriate  selection  of  the  lengths 
and  path  sequence  of  these  strings  enabled  one  to  achieve  the  desired  geo¬ 
metrical  shape  of  the  ensemble  when  it  was  hoisted  into  the  incident  beam; 
both  the  upward  and  downward  pull  of  the  nylon  threads  in  the  orientation 
mechanism  itself  tightened  these  strings  to  hold  the  individual  spheres  in 
place.  Since  each  sphere  could  be  slid  a  little  along  the  strings  passing 
through  its  own  eyelets,  the  precise  relative  positions  of  spheres  were 
adjusted  manually.  This  practice  worked  particularly  well  for  2  spheres 
in  a  discrete  set  of  mutual  separations  and  for  4  to  8  spheres  forming  a 
contacting  square  or  cubic  array,  and  such  arrays  could  be  rotated 
azimuthally  without  deformation. 

Continuous  separation  of  2  spheres  for  the  0=0°  measurements  were  per¬ 
formed  only  for  the  orientation  where  2  spheres  were  aligned  along  the  in¬ 
cident  direction.  In  this  case  the  first  sphere  was  positioned  in  the  normal 
target  site  of  the  orientation  mechanism.  The  second  sphere  was  suspended, 
like  a  pendulum  bob,  from  the  pivoting  point  of  the  orientation  mechanism 
near  the  ceiling  by  means  of  a  separate  nylon  string.  Another  thin,  long, 
nylon  string  was  passed  through  2  eyelets  of  this  second  sphere,  the  string 
ends  separately  going  through  two  floor-level  eyelets  straddling  (but  "down¬ 
stream"  from)  the  incident  beam.  This  string  also  served  as  a  stabilizing 
agent  of  the  movable  sphere,  and  the  continuous  separation  of  2  spheres  was 
accomplished  by  simultaneously  pulling  both  string  ends. 

These  target  separation/orientation  techniques  have  remarkable  simplicity, 
speed,  and  reliability  if  the  following  precautions  are  taken:  (l)  The  en¬ 
semble  should  stay  within  the  acceptable  region  of  the  incident  beam.  Thus , 
the  separation  distance  between  2  spheres  can  not  be  indefinitely  large. 

(2)  For  the  discrete  2-sphere  separation  case,  the  target  height  is  depen¬ 
dent  on  the  separation  distance.  (3)  For  the  continuous  separation  case, 
it  is  somewhat  difficult  to  hold  both  spheres  at  the  same  height.  The  sphere 
being  pulled  was  observed  to  go  higher  than  the  other  by  as  much  as  its  own 
radius  at  a  separation  distance  of  five  sphere-diameters  if  care  was  not 
exercised  during  the  pull.  (A  servo  controlled  motor-gear  mechanism  was 
later  built  which  seemed  to  remedy  this  difficulty.) 

We  have  to  remark  that  such  an  ensemble  of  string-assembled  multiple 
spheres  can  only  be  rotated  azimuthally  without  geometrical  deformation 
using  the  present  orientation  mechanism. 

2.4  More  Remarks  on  Antennas  and  Sources  of  Experimental  Errors 

We  have  already  reported  the  important  role  of  antennas  in  indoor 
precision  measurements  (Wang  et  al.,  1977;  Wang  and  Greenberg,  1978).  In 
combination  with  a  good  anechoic  chamber,  the  antenna’s  design  not  only 
determines  the  size  range  of  scatterers  that  can  be  investigated,  but 
can  also  be  used  to  reduce  the  unwanted  background  radiation.  Careful  align¬ 
ment  of  the  transmitting  and  receiving  antennas  insures  the  syrametricalness 
of  the  radiation  and  reception  patterns  about  the  beam  axes.  During  the 
period  when  the  data  in  this  report  was  taken,  the  antenna  alignment  could 
only  be  carried  out  by  adjusting  the  orientations  of  both  transmitting  and 


receiving  apertures  to  yield  the  maximum  received  signal.  Subsequent  fine 
adjustments  were  made  to  see  if  at  0=0°  the  amplitude-phase  plot  of  a  small 
sphere  would  retrace  as  the  sphere  was  moved  across  the  beam.  If  we  changed 
both  the  incident  and  receiving  polarizations  through  90°  after  this  adjust¬ 
ment  was  made,  however,  an  unacceptable  retrace  was  observed.  The  mechanical 
difficulty  of  precisely  keeping  the  dipole-disk  assemblies  in  the  geometrical 
centers  of  each  antenna  as  they  were  rotated  explains  why  both  polarizations 
were  kept  fixed  (vertical)  throughout  this  experiment. 

Antenna-target  multiple  scattering  is  another  source  of  experimental 
error  for  0=0°  measurement  and  causes  systematic  errors  even  in  single¬ 
sphere  runs.  In  addition  to  the  incident  plus  singly  scattered  waves  arriving 
at  the  receiver  antenna,  we  have  2  first-order  multiple  scattered  signals, 
neglecting  the  smaller  higher  order  multiple  reflections.  These  2  signals 
go  through  the  path  sequences  transmitter-scatterer-transmitter-receiver  and 
transmitter-receiver-scatterer- receiver,  respectively.  Analyses  (Lind,  1966; 
Wang,  1968)  showed  that  the  percentage  error  of  the  amplitude  and  phase 
measurements  at  0=0°  increased  linearly  with  the  ratio  |s(ir) |/|s(0) |  o'  the 
target,  decreased  linearly  with  the  antenna-target  separation,  and  was  pro¬ 
portional  to  the  backscatter  cross  sections  of  antennas.  With  the  latter 
two  parameters  fixed,  the  error  is  seen  to  be  less  for  targets  with  smaller 
|s(ir) | / | S( 0 ) |ratio  like  the  softer  particles  made  from  expanded  polystyrene 
and  larger  for  harder  particles  prepared  from  acrylic  material.  Errors  in 
the  latter  case  may  be  as  large  as  ''*1055  in  amplitude  and  M.20  in  phase).  A 
general  remedy  for  this  source  of  error  is  not  yet  well  established. 

A  controversial  requirement  on  the  minimum  antenna- target  separation 
distance  (Silver,  19**9,  1962;  Beard  et  al.  ,  1962;  Hansen  and  Bailin,  1959; 
Rhodes,  195*0  did  not  seen  to  have  a  critical  effect  on  the  measurements. 

We  employed  only  about  half  of  the  far-zone  distance  in  this  investigation. 

3.  EXPERIMENTAL  RESULTS 
3.1  Forward  Scattering  (8=0°) 

Each  P,  Q  plot  in  this  section  is  a  cartesian  representation  of  the 
complex  forward  scattering  amplitude  S(0)  as  a  function  of  orientation  angle 
X  or  the  mutual  separation  ks=2irs/A  of  the  particle  ensemble  (see  Fig.  2). 

In  the  complex  plane  the  dimensionless  P  and  Q  components  are: 

P  =  Im  {S(o)>  ,  Q  “  — -*■  Re  {S( 0 ) )  ,  (l) 

k2G  k2G 


where  G  is  the  appropriate  geometrical  cross  section  of  the  ensemble  (the  sum 
of  geometrical  cross  sections  of  the  component  spheres),  and  Q  is  the  so  called 
"extinction  efficiency".  A  vector  drawn  from  the  coordinate  origin  to  each  x 
(or  ks)  position  along  the  curve  yields  the  complex  value  of  S(0)  at  x  (or  ks). 
The  phase  shift  $(0)  of  the  0=0°  scattered  wave  is  given  by  the  angle  between 
this  vector  and  the  P  axis,  while  the  projection  of  this  vector  into  the  cali¬ 
brated  Q  axis  gives  the  extinction  efficiency  (Cext/G).  The  length  of  this 
S(0)  vector  represents  the  absolute  value  |s(0)|,  and  its  numerical  value  is 
obtained  by  comparing  this  length  with  that  of  the  "standard"  or  calibration 


vector  (obtained  from  a  sphere)  provided  in  each  plot.  The  numerical  value 
of  |s(0) |  for  the  calibration  sphere  is  given  in  Table  1. 

The  P,  Q  plots  in  this  subsection  are  subdivided  according  to  number, 
size,  refractive  index  and  method  of  separation  of  the  component  spheres, 
the  parameters  of  which  are  listed  in  each  P,  Q  plot  and  also  tabulated  in 
Table  1  along  with  relevant  Mie  scattering  quantities. 

3.1.1  Arrays  of  2  Identical  Spheres  in  Discrete  Steps  of  Mutual  Separation 
and  of  Rotation  in  the  Incident  k-H  Plane 

Figures  3A-3I,  1*A-4h,  and  5A-5F  are  the  experimental  P,  Q  plots  for 
3  pairs  of  2-sphere  ensembles.  The  pairs  differ  from  one  another  only  in 
the  size  of  the  component  spheres.  Target  parameters  and  the  mutual  sep¬ 
aration  ks  is  shown  in  each  plot.  At  the  orientation  k,  the  array  is  aligned 
on  the  incident  direction.  It  is  then  continuously  rotated  through  90°  in 
the  k-H  plane  of  the  incident  wave  to  display  the  continuous  curve,  finally 
arriving  at  the  orientation  H  where  the  array  is  perpendicular  to  the  in¬ 
cident  k  vector.  Fiducial  marks  along  the  curve  denote  the  orientation  angle 
X- 

All  component  spheres  have  refractive  indexes  near  m=1.365  and  resemble 
those  of  water  or  ice  in  the  optical  spectrum.  Such  a  2-sphere  ensemble 
displays  the  simplest-looking  P,  Q  plot  among  those  made  by  other  multiple 
spheres,  but  a  detailed  explanation  is  not  yet  available.  However,  some 
striking  phenomenological  features  are  evident:  (a)  The  phase  shift  4»(0) 
of  the  0=0°  scattered  wave  is  invariably  larger  at  the  k-orientation  than 
at  the  H  orientation  for  these  particle  sizes,  (b)  At  the  H  orientation, 

6(0)  is  about  the  same  as  that  of  a  component  sphere,  but  the  magnitude 
|s(0) |  is  about  twice  as  large;  i.e.,  the  two  spheres  scatter  the  incident 
wave  more  or  less  independently,  (c)  As  the  array  is  rotated  from  k  to  H 
orientation,  the  tip  of  S(0)  vector  generates  a  clockwise  arc  or  spiral, 
converging  toward  the  H  orientation  faster  as  x  increases.  At  x=60°, 
the  S(o)  vector  is  already  near  that  at  x=90°,  the  H  orientation.  Crudely 
speaking,  one  sphere  emerges  out  of  the  other's  shadow,  (d)  As  the  mutual 
separation  ks  increases,  the  above  feature  (c)  is  more  pronounced,  resulting 
in  the  shrinking  of  the  overall  size  of  the  P,  Q  plot.  This  effect  is  a  vivid 
representation  of  the  decrease  in  dependent  scattering. 

3.1.2  Arrays  of  2  Dissimilar  Spheres  in  Discrete  Steps  of  Mutual  Separation 
and  of  Rotation  in  the  Incident  k-H  Plane 

P,  Q  plots  of  an  array  of  2  different  sized  spheres  being  rotated  in  the 
k-H  plane  are  shown  in  Figures  6A-6f  for  six  steps  of  separation  ks.  In 
addition  to  a  standard  calibration  sphere,  the  vector  S(0)  of  each  component 
sphere  was  measured  independently.  The  results  are  shown  in  the  same  P,  Q 
plot  as  2  dotted  lines  marked  SPHERE  1  and  SPHERE  2,  respectively.  All  charac 
teristic  features  in  the  previous  section  3.1.1  apply  here  also.  Notice,  in 
particular,  that  at  H  orientation  the  S(0),  vector  of  the  array  is  nearly  equal 
to  the  vector  sum  of  the  individual  spheres. 

P,  Q  plots  for  another  2-sphere  array  whose  component  spheres  are  nearly 
the  same  in  size  but  different  in  refractive  index  are  shown  in  Figures  7A-7H 
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if  8»0°  scattering  "by  an  array  of  2  identical  spheres  as  the 
continuously  varied.  The  separate  graphs  refer  to  discrete 
See  sections  3.1  and  3.1.1. 
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Figs.  6E  and  6f  Same  as  Figs.  6A-6d 
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Same  as  Figs.  6A-6D,  except  the  2  spheres  here  are  dissimilar  in  their  refractive  indexes 


spheres  here  are  dissimilar  in  their  refractive  indexes 
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TABLE  2  COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL  |S(0)|  s  AND  $(0)  s 

The  experimental  8=0°  scattering  amplitude  |s(0) |  and  phase  shift  $(0)  for 
each  contacting  2-sphere  array  at  the  broadside  incidence  (ll)  orientation 
is  compared  with  the  corresponding  theoretical  predictions  resulting  from 
the  summation  of  independent  scattering  by  these  2  Mie  spheres.  All  spheres 
are  of  acrylic  material,  with  refractive  index  m=1.6l0-i0.004. 


Array  x  Experimental  Theoretical 


No. 

Size  of 

Component 

Sphere 

|S(0)| 

*(o) 

|S(0)1 

*(0) 

1 

1.237 

1.78 

21.2° 

1.689 

17.3° 

2 

1.368 

1.93 

22.5° 

2.338 

20.7° 

3 

1.856 

5.40 

33.5° 

6.030 

41.5° 

4 

2.166 

8.49 

46.2° 

8.504 

48.2° 

5 

3.085 

36.25 

53-0° 

20.16 

71.7° 

6 

3.733 

28.15 

93.8° 

27.74 

91.2° 

7 

4.341 

38.68 

106.6° 

35.26 

106.4° 

8 

4.980 

41.75 

115.0° 

39-06 

113. 9° 
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for  8  steps  of  center-to- center  separation  ks.  In  the  optical  spectrum, 
the  refractive  index  of  one  component  sphere  resembles  that  of  ice  or 
water,  while  the  other  that  of  silicates.  Identical  P,  Q  plots  were  ob¬ 
served  in  both  k-H  plane  rotations  whether  the  silicate  particle  was  in¬ 
itially  (at  k  orientation)  placed  in  the  shadow  of  the  ice  particle  or 
vice  versa.  Here  again  the  conspicuous  feature  of  the  minimum  dependent 
scattering  effect  is  obvious  at  the  H  orientation;  the  resultant  vector 
S(0)  falls  in  the  vicinity  of  the  vector  sum  of  the  S(0)  vectors  of  the 
(independent)  component  spheres  at  this  orientation.  Characteristics  (a), 
(c)  and  (d)  described  in  Section  3.1.1  are  also  observed  here. 

3.1.3  Contacting  2  Identical  Spheres  of  Eight  Array  Sizes 

The  scattering  by  8  pairs  of  2  contacting  identical  spheres  of  acrylic 
material,  each  pair  differing  from  the  others  only  in  size,  were  measured 
separately  to  give  the  signatures  shown  in  Figures  8A-8H.  Notice  that  each 
pair  has  the  same  volume  as  that  of  a  prolate  spheroid  of  elongation  2:1 
whose  semi-minor  axis  is  equal  to  the  sphere  radius.  Due  to  the  high  re¬ 
fractive  index  of  these  particles,  corresponding  to  that  of  silicate  par¬ 
ticulates  in  the  optical  spectrum,  we  also  notice  marked  differences  from 
the  P,  Q  plots  of  dylite  (expanded  polystyrene)  particles  in  preceding 
sections.  It  was  rather  difficult  to  prepare  more  than  two  identical 
spheres  from  the  commercially  machined  sphere  stock,  and  it  was  therefore 
impractical  to  do  the  (independent)  component  sphere  runs  quickly  enough 
before  the  microwave  compensation  drift  became  perceptible.  In  addition, 
these  silicate-type  spheres  have  a  high  back  to  forward  scattering  ratio 
(i(n)/i(0))  as  indicated  by  Table  1,  and  hence  the  measurements  may  contain 
pronounced  errors  due  to  the  antenna-target  multiple  scatter  effect  (Lind, 
1966;  Wang,  1968).  We  could  not  correct  this  effect  by  performing  the  then 
purported  array-pulling  technique  without  seriously  deforming  the  array 
geometry.  Hence,  we  reproduce  here  the  direct  observations  without  such 
technical  modifications. 

A  closer  analysis  of  these  P,  Q  plots  and  a  comparison  with  those  by 
ice-like  particles  in  the  preceding  sections  reveal,  however,  the  following 
similarities  and  contrasts:  (a)  At  the  H  orientation,  where  an  array 
presents  its  broadest  face  with  respect  to  the  incident  direction,  S(0)  of 
the  array  is  close  to  that  of  the  vector  addition  of  S(0)  vectors  of  com¬ 
ponent  spheres,  with  the  exception  of  two  particular  arrays  shown  in  Figures 
8c  and  8e.  This  vector  addition  of  independent  scattering  as  obtained  by 
Mie  theory  is  compared  with  the  observations  as  shown  in  Table  2.  (b)  For 

smaller  array  sizes.  Figures  8A-8E,  the  tip  of  the  S(0)  vector  generates  a 
clockwise  arc/spiral  as  the  array  is  rotated  from  k  to  H  orientation.  As 
the  size  increases  further,  this  trend  disappears,  resulting  in  a  complex 
variation  with  respect  to  the  changing  orientation.  Such  a  P,  Q  plot  bears 
little  resemblance  to  that  obtained  from  a  similarly  sized  array  of  ice¬ 
like  particles  with  lower  index  of  refraction,  (c)  For  all  sizes,  both  the 
magnitude  |s(0)|  and  phase  $(0)  change  most  rapidly  during  the  orientation 
change  10®sxs60°,  and  at  x=60°  S(0)  already  converges  to  that  near  x=90°, 
the  H  orientation,  (d)  If  an  array  is  properly  sized  as  in  Figure  8f,  the 
total  cross  section  of  the  array  may  vary  by  a  factor  *v7  during  such  a  k-H 
plane  rotation. 
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Figs.  8A-8D  P,3  plots  of  0=0°  scattering  by  arrays  of  contacting 

2  identical  spheres  as  the  array  orientation  (x)  is  varied.  The 
separate  graphs  refer  to  arrays  of  different  size.  See  sections 
3.1  and  3.1.3. 


3.1.4  Continuous  Separation  of  2  Identical  Spheres  Along  Incident  Direction 

In  this  section,  we  present  a  detailed  investigation  of  the  dependent 
scattering  process  when  2  spheres  are  continuously  pulled  apart  along  the 
direction  of  the  incident  radiation,  thereby  recording  a  P,  Q  plot  as  a 
function  of  separation.  In  this  special  orientation,  one  sphere  is  always 
in  the  geometrical  shadow  of  the  other,  and  a  more  coherent  understanding 
of  their  mutual  dependence  may  be  possible.  This  was  done  for  4  pairs  of 
spheres,  one  of  which  has  ice-like  particles  (Fig.  9A) ,  while  others  re¬ 
semble  silicates  (Figs.  9B-9C).  Target  parameters  can  be  found  in  each 
figure  as  well  as  in  Table  1.  The  running  numbers  in  each  P,  Q  plot  denote 
the  mutual  separation  ks=2irs/X,  where  s  is  the  center-to-center  distance 
between  two  spheres. 

Some  striking  phenomena  observed  in  these  P,  Q  plots  are:  (a)  In 
general,  the  tip  of  the  S(0)  vector  travels  in  a  counterclockwise  curve 
in  each  P,  Q  plot  as  ks  increases,  converging  toward  a  limit  corresponding 
to  ks=®  with  ever  decreasing  speed.  This  limit  represents  the  summation 
of  independent  scattering  by  two  spheres,  as  shown  in  each  plot  by  an  ex¬ 
tended  straight  section,  half  of  which  (the  solid  section)  being  the  S(0) 
vector  of  an  isolated  single  sphere  as  measured  during  the  same  experimental 
run.  (b)  If  the  center-to-center  distance  between  two  spheres  (s)  is  less 
than  about  1.5  diameters  of  a  sphere,  the  P,  Q  vs  s  plot  is  a  featureless 
curve.  Beyond  this,  however,  the  plot  exhibits  wavy  structures  superimposed 
on  a  main  locus,  indicating  the  S(0)  vector  of  the  array  oscillates  about  a 
median  function  of  s.  (c)  The  extinction  cross  section  of  the  array  is  not 
always  a  monotonous  function  of  s.  Depending  on  the  size  and  refractive 
index  of  component  particles,  there  may  be  a  major  turning  point  (in  addition 
to  small  wiggling  variations)  around  which  Cg^r  reverses  its  increasing/de¬ 
creasing  trend  as  the  separation  proceeds,  (d)  Even  more  subtle  is  the 
fact  that  an  array's  Cgyip  does  not  necessarily  reach  its  maximum  or  minimum 
value  when  the  spheres  are  in  contact  or  are  well  separated.  The  minimum/ 
maximum  value  of  CEXT  is  highly  dependent  on  the  size  and  refractive  index 
of  the  component  particles  as  well  as  on  their  mutual  separation.  In  one 
extreme  case,  Cgx»r  does  not  vary  appreciably  as  a  function  of  ks  (Figs.  3A-3I), 
while  for  another  case,  it  changes  by  a  factor  of  7  (Fig.  9B).  Furthermore, 
no  simple  relationships  seem  to  exist  between  maximum  or  minimum  Cext'e  an<* 
component  sphere's  Cgy^* 

3.1.5  Multiple  Spheres  in  Contact 

These  multiple  spheres  were  prepared  by  molding  expandable  polystyrene 
in  three  different  size  stainless  steel  molds,  in  order  to  obtain  at  least 
nine  identical  spheres  for  each  size.  0=0°  measurements  for  all  of  them 
were  carried  out  to  insure  they  possessed  nearly  the  same  |s(0)|  and  $(0) 
in  each  size  group.  Four  or  eight  spheres  in  each  group  were  assembled  on 
nylon  strings  to  form  arrays  of  contacting  square  or  cubic  geometry,  and 
the  arrays  were  rotated  in  the  incident  k-H  plane  to  obtain  P,  Q  plots  such 
as  those  shown  in  Figures  10A-10C  or  in  Figures  11A-11C.  The  size  parameters 
(x's)  of  the  component  spheres  are  x=3.120,  3.752  and  4.678,  respectively, 
and  the  corresponding  Mie  scattering  quantities  are  tabulated  in  Table  1. 
Numerical  data  of  the  measured  amplitude  |s(0)|,  the  phase  shift  $(0)  and 
extinction  efficiency  Qgjnp  are  listed  in  Table  3. 


Fig.  9A  P,Q  plot  of  9*0°  scattering  by  an  array  of  2  identical  spheres  in  continuous 
separation  along  the  incident  direction.  The  curve  is  parameterized  by  the  mutual 
separation  (ks).  See  sections  3.1  and  3.1.^. 


TABLE  3  FORWARD  SCATTERING  QUANTITIES  OF  CONTACTING  MULTIPLE  SPHERES 


The  orientation  dependence  of  the  0=0°  scattering  quantities  |s(0)|,  $(0)  and  QEXp  for 
arrays  of  2n  (n=l,2,3)  contacting  identical  spheres  are  listed  below  for  3  different 
size  groups.  In  each  array,  all  component  spheres  are  assembled  on  nylon  strings  to 
form  a  dumbbell  (n»l),  a  square  (n=2)  or  a  cube  (n=3),  respectively.  The  array  is 
then  continuously  rotated  in  the  incident-wave  k-H  plane.  At  the  orientation  x=0° 

£  straight  line  passing  through  the  centers  of  adjoining  two  spheres  is  parallel  to 
k  but  perpendicular  to  H,  while  at  X“90°  it  is  perpendicular  to  k  but  parallel  to  H. 
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♦Only  for  this  array,  the  two  spheres  were  in  slight  separation.  The  observed  gap 
distance  was  0.15  cm.  which  gives  ks=9-65  • 


Figs.  10A — OC  Three  P,Q  plots,  each  of  which  shows  the 

orientation  (x)  dependence  of  S(0)  for  a  contacting  square 
array  of  4  identical  spheres.  At  xa0°  (or  90°),  one  side  of 
the  square  is  parallel  to  the  incident  k  vector,  and  the 
square  is  then  continuously  rotated  in  the  k-H  plane  to 
generate  the  plot.  See  sections  3.1  and  3.1.5. 


We  have  previously  reported  a  few  conspicuous  extinction  signatures 
exhibited  by  such  ensembles  (Wang,  1980).  Additional  findings  include: 

(a)  An  array  of  regular  geometry  (suggestive  of  atoms  in  a  crystal  lattice) 
shows  a  large  variation  of  |s(0)|  and  $(0)  with  respect  to  its  orientation 
change.  Although  not  shown  in  this  report,  a  randomly  assembled  array  con¬ 
sisting  of  the  same  component  spheres  yields  much  smaller  variations  of 
| S ( 0 ) |  and  $(0).  (b)  As  expected  the  P,  Q  plots  appear  almost  symmetric 

about  x=^5°  with  respect  to  changes  in  x*  The  apparent  deviation  from  this 
perfect  symmetry,  although  slight,  is  presently  attributed  to  a  slight 
variation  in  target  parameters  from  one  component  sphere  to  another,  along 
with  a  slight  distortion  of  array  geometry  arising  from  the  manual  align¬ 
ment  procedure.  This  indicates  the  critical  dependence  of  |s(0)j  and  $(0) 
on  these  parameters  even  for  0=0°  measurements,  (c)  Note  the  comparisons 
between  Figures  10A  and  11A,  10B  and  11B,  and  10C  and  11C.  Each  pair 
represents  two  P,  Q  plots  of  4  and  8  identical  spheres,  respectively,  differing 
only  on  the  total  number  of  particles  in  the  array.  The  signatures  are 
similar  except  that  |s(0)|  for  the  8-sphere  ensemble  is  about  twice  the 
value  of  the  4-sphere  ensemble.  Detailed  numerical  data  are  also  tabulated 
in  Table  3.  This  implies  no  appreciable  interference  between  two  4-sphere 
arrays  in  an  8-sphere  ensemble  that  are  parallel  to  the  k-H  plane,  and 
hence  these  two  subarrays  are  scattering  independently.  A  generalization 
of  finding  (c)  may  be  inferred  from  comparisons  in  Table  3.  (d)  At  the 

particular  orientation  x=0°»  "the  2-,  4-  or  8-sphere  ensembles  in  each  size 
group  have  practically  the  same  phase  shift  $(0)  but  |s(0)|  is  nearly  pro¬ 
portional  to  the  total  number  of  component  spheres.  One  is  tempted  to  say 
these  aggregated  spheres  can  be  split  in  subgroups ,  all  parallel  to  the  in¬ 
cident  direction,  each  scattering  independently  of  others. 

3.2  Side  Scattering  and  Angular  Distribution 

Complementary  to  extinction  measurements ,  observation  of  scattered 
waves  from  scattering  angles  other  than  6=0°  provides  an  additional  wealth 
of  information  about  the  scattering  process.  Because  of  the  difficulties 
experienced  in  performing  angular  experiments  (Wang  and  Greenberg,  1978), 
the  0=90°  scattering  measurement  was  first  tried  before  the  more  extensive 
angular  distribution  studies.  At  0=90°,  the  compensation  of  unwanted  back¬ 
ground  radiation  can  be  neglected  due  to  its  smallness  in  comparison  with 
the  true  scattered  signal,  and  the  target-to-wall-to-antenna  secondary 
scattering  could  be  reduced  by  placing  good  absorbers  on  that  portion  of 
wall  facing  the  receiver  antenna.  The  resulting  measurements  reveal 
simple-looking  yet  specular  intensity  variation  as  a  pair  of  spheres  was 
rotated  in  the  beam.  Encouraged  by  these  results ,  ve  measured  the  angular 
distribution  of  i 1 ( 9 )  over  the  range  4O°^0il4O°  in  which  the  background 
level  was  relatively  low. 

3.2.1  Side  Scattering  (0a9O°) 

Figures  12A  and  12B  show  the  observed  intensity  variation  as  a  function 
of  ensemble  orientation  x  for  a  pair  of  contacting  identical  spheres  made 
of  expanded  polystyrene.  As  x  is  continuously  varied  for  the  display  of 
il(9),  we  notice  a  fairly  symmetrical  intensity  profile  about  x  =  8/2,  which 
confirms  the  symmetry  relation,  Eq.(l*0  of  Section  2.2.  Disregarding  the 
low  signal-level  portion,  i j ( 0 )  goes  through  at  least  three  maxima  and 


minima  vs  x  change ,  presumably  an  interference  phenomenon  between  the 
dependently  and  independently  scattered  waves  from  the  component  spheres. 
Independently  scattered  waves  alone  could  not  account  for  the  precise  x 
positions  where  these  maxima/minima  occurred.  The  most  striking  fact  is 
its  magnitude  at  x=®/2.  In  comparison  with  the  intensity  of  a  single 
sphere,  as  shown  in  each  figure  and  in  Table  1,  this  array  scatters  nearly 
28  times  that  of  a  component  sphere  at  0=90°  and  22  times  at  0=91°.  While 
the  summation  of  independent  scattering  correctly  predicts  the  occurrence 
of  a  maximum  at  this  particular  xs®/2  because  of  the  in-phase  arrivals  of 
scattered  signals  from  these  two  spheres,  the  maximum  intensity  can,  at 
most,  be  4  times  that  of  a  single  sphere  by  this  token. 

Figures  12C-12H  show  the  ij(90)  vs  x  plots  for  six  pairs  of  2-sphere 
ensembles,  in  the  order  of  increasing  particle  size.  Each  pair  consists 
of  2  contacting  identical  spheres  of  acrylic  material  with  index  of  re¬ 
fraction  ra=1.6l0-i0.004.  Each  pair  resembles  a  silicate  dumbb  11  particle 
in  the  optical  region.  This  series  of  measurements  were  performed  to  in¬ 
vestigate  in  more  detail  the  size  dependence  of  specular  scattering  at 
0=90°.  A  few  conspicuous  features  from  these  plots  are:  (a)  Similar  to 
the  previous  pair  of  polystyrene  spheres ,  the  ij  (90)  vs  x  profile  is  symme¬ 
trical  about  x=45°  at  which  the  intensity  is  also  at  a  maximum,  although 
not  necessarily  the  maximum  over  the  entire  x  range  (see  Fig.  12E). 

(b)  This  maximum  is  separated  from  the  adjacent  minimum  (or  maximum)  by 
ever  decreasing  x  intervals  as  the  particle  grows  in  size  -  three  maxima 
and  two  minima  within  45°  rotation  from  x=^5°  for  a  particular  ensemble 
in  Figure  12H.  (c)  In  comparison  with  *1 (90)SNGLgpRERE  of  a  component 

sphere  shown  in  each  plot  and  in  Table  1,  the  onset  of  specular  scattering 
at  which  i)(90)  of  the  array  is  greater  than  4i j (90)siJGLSPK3RE  seems  to 
occur  only  for  moderately  larger  particles  (Figs.  12G  and  12H)  and  ex¬ 
clusively  at  x~1|50.  For  the  particular  ensemble  shown  in  Figure  12H, 
ii(90)=44i1(90)SNGI£pHERE  ! 

3.2.2  Angular  Distribution  (4o°<0<l4o°) 

Angular  distribution  data  for  an  array  of  two  contacting  identical 
spheres  made  of  expanded  polystyrene  is  displayed  in  Figure  13A  for  two 
principal  orientations  k  and  H  in  which  the  array  axis  is  parallel  to  the 
incident  k  and  H  vectors,  respectively.  The  size  of  the  errors  associated 
with  the  smallness  of  the  scattered  intensity  (except  for  6<50°) ,  which  is 
possibly  mixed  with  uncompensated  background  radiation,  with  the  uneveness 
of  the  floor  (which  defines  the  scattering  plane),  and  with  the  mechanical 
alignment  of  the  receiving  antennas  (Wang  and  Greenberg,  1978)  are  difficult 
to  estimate;  we  present  the  data  as  observed. 

The  result  for  another  2-sphere  ensemble  is  shown  in  Figure  13B.  The 
two  identical  spheres  are  of  acrylic  material  and  are  the  same  pair  which 
showed  specular  reflection  at  0=90°  scattering.  More  regularly  shaped 
scattering  lobes  are  observed  for  this  particular  ensemble,  both  for  the  k 
and  H  orientations.  When  the  array  axis  makes  an  orientation  angle  x=**0° 
from  the  incident  direction  toward  the  receiving  direction,  we  see  a  broad 
scattering  lobe  centered  around  0=80°  peaking  at  a  large  intensity.  This 
is  another  confirmation  of  specular  scattering  at  x=®/2  for  this  particular 
array.  Note,  however,  that  the  standard  calibration  sphere  used  in  this 
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Fig.  13A  Angular  distribution  ij(e)  for  an  array  of  contacting 
2  identical  spheres  at  two  special  array  orientations  k  (array 
axis  parallel  to  k)  and  H  (array  axis  parallel  to  H). 
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run  was  later  found  to  be  questionable,  and  hence  the  absolute  magnitude 
calibrated  in  this  Figure  13C  may  be  in  error. 

U.  SUMMARY  REMARKS 

We  have  reported  experimental  results  on  scattering  by  aggregated 
spheres  of  simple  geometrical  shapes.  Notable  findings  through  these 
microwave  measurements  are:  (a)  An  ensemble  of  spheres  of  fixed  geometry 
possesses  intricate  variations  in  its  8*0°  scattering  with  respect  to  the 
orientation  change  of  the  ensemble,  more  regularly  arrayed  ensembles  ex¬ 
hibiting  more  pronounced  variation  and  symmetry,  (b)  When  the  axis  of 
two  spheres  is  perpendicular  to  the  incident  direction,  the  0=0°  scattering 
by  the  ensemble  is,  to  a  good  approximation,  a  simple  summation  of  in¬ 
dependently  scattering  two  spheres.  To  a  lesser  degree,  this  finding 
applies  to  two  arrays  aligned  orthogonally  to  the  incident  direction. 

(c)  The  dependent  scattering  effect  marks  its  presence  at  0=0°  scattering 
if  one  of  the  two  neighboring  spheres  lies  within  the  major  forward  scat¬ 
tering  lobes  of  the  other.  The  angular  extent  of  this  effect  decreases 
nonmonotonically  with  respect  to  the  increasing  mutual  separation,  (d)  The 
P,  Q  plot  of  a  contacting  2-sphere  ensemble  displays  complex  but  systematic 
changes  as  the  aspect  angle  of  the  array  varies.  The  pattern  is  more  com¬ 
plex  for  higher  refracting  particles  than  for  particulates  with  lower  in¬ 
dexes  of  refraction.  A  similar  trend  exists  as  the  size  of  the  component 
spheres  increases.  In  particular,  an  interesting  correlation  seems  to 
exist  between  an  ensemble  of  2  identical  spheres  and  a  prolate  spheroid 
of  2:1  elongation  possessing  equal  volume  and  refractive  index,  (e)  The 
extinction  of  light  by  a  pair  of  spheres  aligned  along  the  incident  direction 
depends  strongly  on  the  size,  index  of  refraction  and  mutual  separation  of 
component  spheres.  It  may  vary  up  to  a  factor  of  7  as  the  separation  in¬ 
creases,  with  ever  slower  convergence  to  the  limit  represented  by  the  sum¬ 
mation  of  extinction  by  component  spheres.  It  is  estimated  that  such  a 
pair  can  be  considered  as  independent  scatterers  if  the  separation  exceeds 
about  10  sphere  diameters,  (f)  Estimates  of  the  extinction  averaged  over  a 
random  orientation  of  an  ensemble  of  multiple  spheres  in  this  size  range 
indicates  that  the  ensemble  obscures  light  more  efficiently  than  a  smooth 
sphere  of  equal  total  volume  (Wang,  1980).  (g)  Certain  contacting  2-sphere 
arrays  in  the  resonance-size  region  scatter  light  specularly  toward  the 
direction  £  when  the  array  axis  bisects  this  direction  and  the  incident 
direction  k. 

Temporal  and  budgetary  problems  did  not  allow  us  to  pursue  another 
aspect  of  important  measurements :  repeating  the  measurements  with  incident- 
wave  polarization  changed  from  vertical  to  horizontal.  Nevertheless,  we 
believe  this  report  contains  very  rare  facts  on  a  number  of  detailed  phen¬ 
omena  about  the  dependent  scattering  of  multiple  spheres. 
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Statement  of  the  problem  studied  -  Many  of  our  microwave  analog  measurements 
of  (mainly)  extinction  and  (some)  angular  scattering  by  irregular  particles, 
taken  over  the  past  two  decades,  have  not  been  published.  Under  this  contract, 
we  have  selected  and  written  up  for  publication  over  sixty  data  sets  which  have 
direct  relevance  to  the  Army's  SMOKE  program  and  which  can  be  obtained  only  by 
the  microwave  analog  method.  These  measurements  provide  experimental  evidence 
on 

(1)  the  role  of  particle  shape  and  particle  inhomogeneity  on  the  ex¬ 
tinction  caused  by  particles  whose  sizes  are  of  the  order  of  the 
incident  wavelength  and 

(2)  the  effect  of  dependent  scattering  by  spheres  which  are  close  to¬ 
gether,  on  both  the  extinction  and  angular  scattering. 

Ifciderstanding  the  physics  involved  in  both  problems  (l)  and  (2)  is  a  major  goal 
of  the  SMOKE  program.  Given  the  solution  of  problem  (l),  one  can  optimize  the 
extinction  of  an  obscuring  cloud  by  a  Judicious  choice  of  the  shape  and  chemical 
make-up  of  its  constituent  particles.  Problem  (2)  addresses  a  different  but 
equally  important  aspect  of  dense  clouds  -  how  far  apart  must  the  individual 
particles  be  so  as  to  treat  them  as  independent  scatterers  of  light?  It  is 
important  here  to  distinguish  between  two  concepts  which  have  a  confusing  and 
overlapping  terminology.  We  do  not  address  in  this  report  "multiple  scattering" 
resulting  from  the  fact  that  a  cloud  is  optically  thick.  (The  physics  of  mul¬ 
tiple  scattering  in  this  sense  is  well  defined  by  the  equation  of  radiative 
transfer).  Rather,  we  address  the  physics  of  dependent/independent  scattering 
(also  referred  to  by  some  authors  as  multiple  scattering)  caused  by  the  close 
proximity  (i.e.  ,  inside  the  "wave  zone")  of  two  or  more  particles.  The  physics 
of  dependent  scattering  is  the  least  understood  and  yet  can  be  the  dominant 
radiative  factor  in  very  dense  clouds. 

Summary  of  the  most  important  results 

1.  Extinction  by  Rough  Spheres  (Represented  by  Stacked  7-Cylinders). 


While  it  is  known  that  particle  shape  dramatically  influences  the  light 
scattering  process  for  large  scattering  angles  (8>90°),  the  effect  of  particle 
shape  on  extinction  is  commonly  held  to  be  non- consequential.  However,  the  re¬ 
sults  of  extinction  measurements  on  23  stacked  7-cylinders  (rough  spheres), 
averaged  over  orientation,  differ  from  the  extinction  produced  by  equal  volume 
spheres  by  up  to  73%  (see  Figs.  7A  and  7B,  p.  49)  in  some  cases. 

2.  Extinction  by  2n  (n=l,2,3)  Multiple  Spheres  as  a 
Function  of  Their  Mutual  Separation. 


When  the 
tinction  (0=0 
scattering  by 


axis  of  2-6pheres  is  perpendicular  to  the  incident  beam,  the  ex¬ 
scattering)  is,  to  a  good  approximation,  a  simple  summation  of 
two  independent  spheres.  To  a  lesser  degree,  this  finding  also 
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applies  to  arrays  aligned  orthogonally  to  the  incident  team.  The  extinction  of 
light  by  a  pair  of  spheres  aligned  along  the  incident  beam  depends  strongly  on 
the  size,  index  of  refraction,  and  mutual  separation  of  the  component  spheres. 
It  may  vary  up  to  a  factor  of  7  as  the  separation  increases ,  with  ever  slower 
convergence  to  the  limit  represented  by  the  summation  of  the  extinction  of 
the  component  spheres.  It  is  estimated  that  such  a  pair  can  be  considered 
independent  scatterers  if  their  separation  exceeds  about  10  sphere  diameters. 

As  a  rule  of  thumb  t  dependent  scattering  becomes  important  for  extinction  when 
the  "back"  particle  is  within  60°  of  the  other's  "shadow"  and  closer  than  10 
sphere  diameters. 

3.  Resonances  in  the  Side  Scattering  of*  Multiple  Spheres. 

When  two  spheres  are  viewed  as  per  the  diagram 
to  the  right,  one  might  anticipate  a  specular  re¬ 
flection.  But  no  theoretical  calculation  had  pre¬ 
dicted  or  experiment  measured  the  extreme  resonances 
described  on  pages  3^-39  of  this  report.  In  one  case 
we  observed  that  the  increase  in  scattered  radiation 
was  a  factor  of  UU  (see  Fig.  12H,  p.  36)  over  that 
produced  by  a  single  sphere. 
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